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ABSTRACT:  In  this  report , a theoretical  study  of  the  radiation 

characteristics  of  rectangular  surface  sources  has  been  made 
when  these  sources  radiate  uniformly  over  their  surfaces  and 
obey  Lambert *s  Cosine  Lav.  This  study  includes  a presentation 
of  the  history  of  the  development  of  such  formulae  up  to  the 
present.  Very  general  equations  of  two  forms  have  been  derived 
giving  the  total  flux  on  an  elementary  receiving  areaT/  when 
this  elementary  receiving  area  has  arbitrary  coordinates  and 
a surface  normal  with  arbitrary  direction  cosines.  The 
classical  method  of  surface  integration  introduced  by  Lambert 
has  bsen  used.  These  equations  are .very  general  in  form,  and 
it  is  shown  how  each  of  the  equations  found  in  the  literature 
become  special  cases  of  these  more  general  equations  in  which 
the  coordinates  of  the  elementary  receiving  area  or  the 
direction  cosines  of  its  surface  normal  have  particular  values. 
A simple  translation  and  rotation  of  coordinates  makes  it 
possible  to  also  consider  the  equally  important  problem  in 
which  the  source  is  permitted  to  have  arbitrary  coordinates 
and  a surface  normal  with  arbitrary  direction  cosines.  Some 
numerical  examples  of  the  radiation  field  of  a square  surface 
source  are  given  and  possible  military  applications  are  cited. 
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RADIATION  CHARACTERISTICS  OP  RECTANGULAR 


SURFACE  SOURCES 


INTRODUCTION 


lo  Often  it  would  tie  quite  useful  If  It  were  possible  to  predict 
the  total  flux  from  a source  of  radiation  of  arbitrary  geometry 
falling  on  a receiver  when  the  source  and  receiver  ara  at  arbi- 
trary distances  from  each  other  and  have  arbitrary  orientations. 

This  flux  might  be  of  any  form  such  as  electromagnetic  flux, 
acoustic  flux,  or  It  may  consist  of  a flow  of  particles.  This 
problem.  In  addition  to  having  many  everyday  applications,  also 
is  one  of  the  problems  facing  the  military  when  devices  are 
constructed  to  detect  potential  targets,  for  certainly  one  of  the 
most  important  problems  in  such  considerations  which  must  be  faced 
is  that  of  determining  the  total  flux  from  the  target  falling  on 
the  receiver o The  receiver,  of  course,  must  also  convert  this 
flux  into  b signs!  having  a sufficiently  high  level  above  the  noise, 
but  such  considerations  form  another  chapter  which  lies  beyond 
the  scope  of  the  ore sent  rbporto  We  shall  be  concerned  here  only 
with  methods  for  determining  the  to^al  flux  falling  on  the  receiver 
and  will  not  be  concerned  with  the  manner  In  which  the  receiver 
utilizes  this  flux. 


Zb  Calculations  of  the  flux  from  sources  of  arbitrary  geometry 
have  not  been  attempted  because  they  ars  usually  regarded  as  too 
complex.  To  be  sure,  they  can  be  carried  out  when  the  source  and 
receiving  element  are  separated  by  distances  which  are  lai'ge 
compared  ©Ith  f; he  dimensions  of  the  source  and  receiving  element 
and  tho  absorption  and  scattering  of  the  Intervening  medium  are 
neglected,  for  in  this  case  tbs  inverse  square  law  can  be  applied. 
There  is  then  no  problem  In  determining  the  total  flux  falling  on 
the  receiver  for  any  position  and  orientation  If  it  Is  known  for 


one  position  and  orientation..  Frequently  however,  receivers  are 


called  upon  to  detect  sources  at  close  range  so  that  the  geometries 
of  both  the  source  and  the  receiver  influence  very  markedly  the 
total  flux  falling  on  the  receiver.  Under  these  conditions* the 
Inverse  square  law  fails,  and  there  is  then  no  simple  way  of 
calculating  the  total  flux  falling  on  the  receiver.,  The  practical 
solution  to  such  a problem  Is  usually  the  engineering  one  In  which 
a receiver  Is  chosen  and  measurements  made0  Such  measurements, 
however,  are  useful  only  under  the  conditions  under  which  they  are 
taken  and  become  of  very  little  value  when  changes  are  made  or  new 
variables  appear. 
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3 o To  predict  the  total  flux  falling  on  a receiver  of  given 
geometry  for  any  position  and  orientation  in  space,  one  must 
first  make  detailed  studies  of  the  radiation  fields  o £ the 
sources  being  considered.  A radiation  field  of  a source  of 
flux  of  some  form  can  be  described  in  various  ways,  usually 
depending  upon  how  it  is  to  be  detected.  Since  we  will  be 
interested  here  in  calculating  the  total flux  falling  on  a 
small  plane  receiving  element  (thermocouple,  photoconducting 
cell,  etc.),  we  will  consider  the  radiation  field  of  such  a 
source  to  be  defined  by  the  total  flux  falling  or.  a small 
elementary  receiving  element  (small  compared  to  other  dimensions) 
placed  at  different  points  in  the  surrounding  region  and  having 
different  orientations « To  determine  the  forms  of  the  radiation 
fields  of  various  sources,  surface  and  volume  integrations 
become  necessary.  Many  of  these  Integrations  have  already  been 
carried  out  for  sources  of  simple  geometries,  providing  certain 
simplifying  assumptions  are  made  concerning  the  properties  of 
the  source  and  flux,  and  the  elementary  receiving  area  is 
placed  in  some  particular  position  or  plane  and  is  oriented 
in  some  particular  direction.  Similar  problems  involving 
"action  at  a distance"  have  also  appeared  in  other  fields  of 
physics  in  which  surface  and  volume  integrations  have  been 
necessary. 


4.  In  the  case  of  surface  sources,  which  will  concern  us  in 
this  report,  some  assumptions  are  usually  made  concerning  the 
nature  of  the  flux  emitted  from  the  surface  in  question. 

Surface  integrations  of  the  tyoe  previously  mentioned  mean 
that  the  contributions  from  each  of  the  surface  elements  are 
simply  added.  Jn  the  case  of  radiant  energy,  if  the  contributions 
of  these  surface  elements  are  to  be  added,  then  the  wave  lengths 
of  the  flux  must  be  both  random  in  phase  across  the  surface 
and  small  compared  to  the  size  of  the  surface  elements  involved. 
This  means  then  that  the  problem  of  evaluating  the  radiation 
field  of  such  a surface  source  emitting  radiant  energy  Is 
reduced  to  one  in  geometrical  optics  rather  than  in  physical 
optics.  Similarly  in  the  case  cf  particles,  the  particles  must 
be  small  compared  to  these  surface  elements  and  they  must  not 
Interact  with  each  other  if  these  surface  Integrations  are  to 
be  of  value.  Assuming  that  these  conditions  are  fulfilled, 
two  additional  assumptions  are  usually  made  concerning  the  nature 
of  the  radiating  surface.  The  first  is  that  the  surface  radiates 
uniformly  over  its  surface  and  the  second  that  it  obeys  Lambert’s 
Cosine  Law  (reference  (a)).  This  of  course  represents  a further 
idealization  of  the  problem,  and  such  a source  is  commonly  referred 
to  as  a "perfect  diffuse  emitter".  Many  sources  will  be  found 
which  do  not  conform  to  such  an  idealization  (reference  (b)). 

Even  though  precautions  are  token  to  be  sure  that  they  radiate 
uniformly  over  their  surfaces,  deviations  still  are  found  because 
Lambert's  Cosine  Law  Is  not  satisfied*  Heated  polished  surfaces 
for  example  always  radiate  more  flux  at  oblique  emergence  than 
would  be  expected  from  Lambert's  Cosine  Law.  Occasionally,  a 
cosine  squared  dependence  with  angle  of  emergence  is  found  for 


2 
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particular  materials  with  particular  surface  conditions o 
Fortunately,  there  are  many  sources  which  approximately  obey 
Lambert’s  Cosine  Law  (e.g.  magnesium  oxide,  finely  ground  glass 
surfaces,  thermal  radiation  from  molten  silver  and  unpolished 
platinum),  particularly  if  the  surface  is  somewhat  rough,  and 
the  range  in  angle  of  emergence  is  not  too  large.  V.'e  will 
therefore  assume  that  the  surface  sources  being  considered  here 
radiate  uniformly  over  their  surface  and  obey  Lambert's  Cosine  Law. 

5.  In  this  report,  we  will  confine  our  attention  to  the  rectan- 
gular emitting  surface.  The  primary  reason  for  this  is  that-  many 
sources  are  rectangular  in  shape.  However,  there  are  also 
sources  which  do  not  possess  simple  forms  and  perhaps  do  not 
radiate  uniformly  over  their  surfaces  or  obey  Lambert' 3 Cosine 
Law.  In  these  cases,  the  source  can  usually  be  divided  up  into 

a series  of  smaller  sources  which  do  approximately  fulfill  these 
conditions  ox  uniformity  and  satisfy  Lambert's  Cosine  Law. 

The  rectangular  source  element  is  one  of  the  most  useful  forms 
for  such  calculations,  and  the  total  flux  falling  on  the  elementary 
receiving  area  from  such  a complex  source  will  simply  be  the  sum 
of  the  individual  contributions  of  each  of  these  smaller  rectan- 
gular sources.  The  general  problem  which  will  concern  us  here 
is  then  more  specifically  that  of  calculating  the  total  flux 
from  the  rectangular  surface  source  (assuming  the  rectangular  ■ 
source  radiates  uniformly  over  its  surface  and  obeys  Lambert's 
Cosine  Law)  falling  on  the  elementary  receiving  area  when  it 
has  arbitrary  coordinates  and  a surface  normal  with  arbitrary 
direction  cosines. 

II.  History 

6.  If  one  assumes  that  the  rectangular  source  radiates  uniformly 
over  its  surface  and  obeys  Lambert's  Cosine  Law,  then  the 
rectangular  emitting  surface  is,  of  course,  a rather  old  problem 
(reference  (c))«  The  origin  of  calculations  on  the  field 
characteristics  of  emitting  rectangular  sources  seems  to  b© 

in  the  papers  of  E.  ?,  Hyue  and  K.  Korden  (reference  (d)) 
published  in  1907  and  1908.  Both  of  these  authors  were  concerned 
with  calculating  the  flux  density  on  an  elementary  receiving 
area  from  a mercury  arc  of  tubular  form,  which  could  be 
approximated  by  a cylinder  or  an  Infinitely  long  strip. 

Apparently  the  first  calculations  for  rectangular  surfaces 
of  finite  dimensions  were  those  made  by  Bordonl  In  Italy  in 
1908.  The  two  cases  in  which  he  was  interested  are  shown 
schematically  In  Fig.  1 and  Fig.  2.  Bordonl  calculated  the 
flux  density  on  an  elementary  receiving  area  located  on  the 
surface  normal  through  the  midpoint  of  one  side  of  the  rectangle 
when  the  receiving  area  was  parallel  to  the  rectangular  surface. 

As  Illustrated  by  Fig.  2,  he  also  concerned  himself  with  the 
case  in \hien  the  surface  normal  to  the  elementary  receiving 
area  was  rotated  until  It  was  parallel  to  one  pair  of  the  sides 
of  the  rectangle.  Both  of  these  equations  are  correct  and  in 
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agreement  with  'the  often  repeated  calculations  for  these  problems 
appearing  In  later  years. 

?•  Another  cf  the  early  workers  on  the  radiation  properties  of 
rectangular  sources  was  B.  Jones  {reference  CD)  In  the  United 
State's  who  published  three  papers  during  the  period  fro*  1909-1911® 
"these  papers  described  sorae  additional  radiation  character  is  tics 
of  rectangular  sources  as  veil  as  other  surface  sources.  The 
particular  problems  involving  rectangular  sources  in  which  he 
was  interested  a, re  those  shown  in  Figs.  3*  4,  5f  and  6.  Figs®  3 
end  4 will  he  reccgni  SO'Cl.  «.s  modifications'  cf  one  of  the  problems 
considered  earlier  by  Bordoni.  The.'  Jones  formula  for  the  flux 
density  for  the  case  represented  by  Fig.  3 was  written  down 
incorrectly  In.  his  first  paper  hut  was  corrected  .In  the'  following 
paper  in  1910*  The  formula  for  the  case  represented  by  Fig.  4 
appears  correctly  in  the  paper  of  1910®  The  problems  suggested 
'by'  Fig.#®  5 and  6 : epresent  the  first  problems  in  which  the 
elementary  receiving  area  was  located  on  surface'  normals  at 
arbitrary  distances  from  the  corner  car  center  of  the  rectangle- 
and  was  oriented  in  directions  other  than  parallel  car  perpen- 
dicular' to  its  edges.  In  the  1909  paper,  the  formula  far  the 
case  represented  by  Fig,  5 was  given  incorrectly.  The  problem 
.suggested  by  Fig.  6 was  solved  correctly  in  Junes*  1911  paper.. 
Unfortunately  Jones  has  made  several  mistakes  so  that  his  papers 
have  been  severely  criticized  by  Is  ter  writers®  These  errors 
on  the  rectangular  source  will,  be  discussed  in  detail  later. 

She  Interest  in  the  rectangular  source  seemed  then  to  subside 
until  1923-  At  that  time  f,  J.  Waldrac  (reference  (g>)  in 
England  published  an  approximate  formula  for  the  flux  density 
on  a horlsontal  surface  when,  the  rectangular  source  was  tilted 
at  an  arbitrary  angle  with  respect  to  the  horizontal.  The 
particular  problem  in  which  he  was  interested  is  shown  in  Fig.  7®. 
to  assumed  that  the-  rectangle'  could,  without  appreciable  error, 
be  placed  on  the  surface  of  a sphere  end  then  considered  the 
flux  density  on  a horizontal  receiving  elsssggt;  pises?!  st 

the  center  of  the  sphere*  This  f omuls  is  then  valid  only  when 
the  distance  fro®  the  rectangular  source  is  large  compared  with 
its  dimensions • 


8.  The  next  real  advance-  in  predicting  the  radiation',  fields  of 
rectangular  sources  was  made  in  Japan  by  25®  Ta.saa.outi  (reference, 
(b))  in  1924®  Yamanouti  calculated  the  flux  density  on  an 
elementary  receiving  area  when  it  had  arbitrary  coordinates 
«o»  ?o,.zo>  and  was  oriented  In  the  three  directions  of  the 
coordinates.  For  these  calculations,  tie  rectangular  source 
was  confined  to  one  of  the  coordinate  planes  at  the  origin  of 
■the-  coordinates  as  suggested  'by  Fig,  8,  The  three  directions 
of  the  surface  normals  of  the  elementary  receiving  area  are 
indicated  in  the  figure.  Taranouti 1 9 considerations  represent 
an.  extensloin  of  the  work  cf  Jones  In  which  two  arbitrary 
ccKjrdinates  were-  possible  for  -one  particular  orientation  of' 
the  receiving  element  {see  Fig®  5).  This  work,  however,  still 


* 
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had  aaore  general  applications,  for  the  three  forzmlae  obtained 
were  looked  upon  as  the  three'  components  of  a "light  vector w. 

To  obtain  an  equation  for  the  flux  density  on  an  elenesr-tery 
receiving  area,  having  arbitrary  coordinates  ana  a surfa.ce  normal 
with  arbitrary  direction  cosines,  Yananccti  suggested  that  a 
scalar  product  could  he  formed  with  the  unit  vector  along  the 
direction  of  the  surface  normal  of  the  elementary  receiving 
area,  but  this  equation  was  never  actually  written  because  of 
the  author* s interest  in  'the  simpler  cases.  The  three  formulae 
obtained  by  Tax.en.outi  are  considerably  more  complex  than  the 
previous  formulae  derived  by  others,  and  some  mistakes  in  signs 
appear  which  will  be  mentioned  later, 

9.  The  next  contribution  to  the  subject  was  that  of  S,  B.  Higble 
and  his  cowarkers  (reference  CD)  in  this  country  in  192 5 and 
1926.  Eighta  gave  a.  formula  for  the  flux  density  an  an  elementary 
receiving  area  with  arbitrary"  coord inaves,  when,  it  was  oriental 
p append icularly  to  the  rectangular  surface.  The  particular 
problem  and  the.  terminology  used  in.  setting  up  this  faonmla  is 
that  'indicated  in  Fig.  9.  Clearly  this  is  equivalent  to  one  of 
the  special  cases  treated  earlier  by  lamanouti.  and  Yaaanouti,s 
corresponding  equation  com  be  reduced  to  Higble *s  equation  when, 
both  are  expressed  in  the'  same  terms,  (assuming  that  the  previously 
mentioned  errors  in  signs  in.  Tasanouti * s equations  have  been, 
corrected).  Higble  also  discussed  two  other  .special  cases  which 
are  of  interest  and  are  shewn  schematically  in  Figs.  10  and  11.. 

The  problem  suggested  by  Fig.  10  is  also  a special  case  of' 
YamancmtiIs  equations  and  can  (assuming  that  the  previously 
mentioned  errors  In  sign  have  been  corrected)  be  brought  into 
agreement  with  YamanoutDs  corresponding  equation  when  both 
equations  are  expressed  in  the  same  terms.  Fig.  11  is  another 
modification  cvr  the  Jones  problem  (see  Fig..  6)  and  Is  & special, 
case  of  the  i*.  '*■  general  equation  in  the  present  paper.  These 
ssss.ti.'Sss  will  ss  dlSEssssst  rors?  Is  dst.sl.X  letter.  Xu  X92® 

0*  Seibert  (reference  (J)}  in  Germany  published  a paper 'in  which 
he  again,  derived,  'the  two  earlier  equations  of  B<orddoi.  In  this 
paper,  however , he  used  spherical  trigonometry  and  projection 
methods  rather  than  the  usual  method  of  surface  integration.  • 
Vilhelm  Kusselt  (reference  (k))  in  Germany  during  the  same  year 
also  published  a paper  showing  how  projection  methods  on  a unit 
sphere  were  equivalent  to  surface  integration. 

10*  P.  Moon  (reference  (1))  and  B,  S.  Spencer  (reference  (m>) 
in  this  country  have  also  carried  out  extensive  calculations  on 
rectangular  sources  using  spherical  trigononetry  and  vector  methods. 
Like  Yamnouti  they  have  been  concerned  with  evaluating  the  flux 
density  on  an.  elementary  receiving  area,,  when  it  has  arbitrary 
■coordinates  and  is  oriented  in  the  directions  of  the  three 
coordinates  (see-  rig.  8).  The  equations  obtained  are  in  slightly 
different  forms-  from  those  of  Yemancuti,  but  can,  after  a.  slight 
faw»fermt5.?m9  be  shown  to  be  equivalent  (except  for  the  errors 
in  sign  in  Xanauouti’s  equations). 
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11*  The  most  recent  papers  on  this  subject  are  those  of 
Cx*  Beths 1 s (reference  Cn>)  In  Germany.  These  papers  present 
a vary  general  equation  for  the  rectangular  source,'  yfaea  the 
source  Itself  Is  confined  to  one  of  the  coordinate  planes  and 
the  elementary  receiving  area  Is  located  on  one  of  the  coordinate 
axes  and  has  arbitrary  orientations*  This  equation  Is  very  long 
and  inquires  over  half  of  an  ordinary  sine  Journal  page  to  record 
It,  'This  equation,  however,  represents  a further  generalisation 
of  the  problem  beyond  the  equations  of  Yamanouti  and  will,  after 
a translation  of  coordinates, 'be  found  to  be  equivalent  to  the 
sore  condensed,  equation  in  the  present  paper  when  the  two  are 
reduced  to  the  same  terminology.  The  author  then  presents  two 
special  cases  of  this  equation  which  are  new  and  will  be  discussed, 
later. 

12.  In  view  of  the  somewhat  unsatisfactory  state  in.  which  -the 
rectangular  surface  source  has  been  left,  it  is  felt  that  it 
sight  be  worthwhile  to  solve  this  problem  independently  to  sea 
If'  the  equations  of  others  can  be'  verified  and  perhaps  add  new 
ana.  simpler  formulae » It  was  also  felt  that  'there  was  a genuine 
need  to  sort  of  coordinate  or  bring  together  all  these  apparently 
very  different  formulae,  so  that  it  would  be  possible  to  see 
wherein,  they  differ,  or  perhaps  are  similar,  and  also  whether 
or  not  they  are  correct,  toother  motivating  force  for  'this 
report  has  been  the  desire  to  calculate  the  radiation  fields 
of  Irregular  surface  sources  which  could  be  approximated  by  a. 
series  of  smeller  irregularly  oriented  rectangular  emitting 
surfaces.  To  do  these  three  things  It  became  clear  that  a very 
general  approach  would  be  preferable  in  which  a very  general 
equation  was  found  which  described  all  these  cases.  He  there- 
fore undertook  to  solve  the  problem  of  calculating  the  total 
flax  from',  a rectangular  emitting  surface  falling  on  an  elementary 
receiving  area  hating  arbitrary  coordinates  aid  a.  surface  normal 
with  arbitrary  direction  cc-slsss.  This  prbbL.es  was  solved  before 
we  found  the  papers  of  Be  the.  It  will  bs  recognised  that  Bethe*s 
problem  is  very  close  to  the  problem,  proposed  here,  far  he  has 
limited  his  rectangular  source  to  one  of  the  coordinate  planes 
and  permitted  his  elementary  receiving  area  to  have  arbitrary 
orientations  along  one  of  the  coordinate  axes,  while  we  are 
suggesting  that  the  rectangular  source  be  limited  to  one  of'  the 
coordinate  planes  at  .the  center  of  the  coordinates  and  the 
elementary  receiving  area  permitted  to  have  arbitrary  coordinates 
and  a surface  normal  with  arbitrary  direction,  cosines.  Thus, 
although  the  integral  equations  for1  the  two  problems  are. 
different,  the  resulting  equations  should  be  reducible  to  each 
other  when  the  two  equations  are  expressed  in  the  same  terms, 
for  in  going  from,  one  problem  to  the  other,  a simple  translation 
of  coordinates,  is  necessary.  This  is  fortunate,  since  each 
of  the  two  sets  of  rather  detailed  calcu3.at.lons  gives  added 
weight  to  the  other • ¥e  should  like  then  in  the  following  to 
give  a rather  brief  derivation  .of  our  formula  for  the  rectan- 
gular source,  when  the  elementary  receiving  area  has  arbitrary 
coordinates  and  a surface  normal,  with  arbitrary  direction  cosines. 


6 


KAVGHD  Report  29®© 


We  will  shear  how  all  the  previously  mentioned  equations  become 
special,  cases  of  this  general  foraaila  and  point  out  werein  we 
think  errors  exist.  Some  examples  will  then  he  giver-  of  radiation 
fields  cf  a square  surface  source.  The  more  complex  calculations 
were  nsde  possible  by  utilising  'the  mechanised  computing 
facilities  of  the  If  ami  Ordnance  Laboratory. 

III.  General  Equations  for  Radiation  Field  of 
Rectangular  Surface  Source 

13*  To  describe  in  core  detail  the  problem  in  which  we  shall 
be  interested,  it  will  be  necessary  to  refer  to  Fig,  12.  In  the 
(X.S)  plane  of  the  coordinate  system,  we  Swire  indicated  schenatf- 
'Sslly  a plane  rectangular  surface  whose  center  is  at  the  origin, 
of  the  coordinate  system.,  tfe  will  assume  that  this  rectangular 
surface  acts  as  a source  of  flux  and  will  not  be  concerned  with 
the  mechanism  by  which  this  flux  arises  • It  is  also  uaiaportant 
as  to  the  form  which  this  flux  takes.  As  suggested  earlier,  it 
may  be  either  a wave  disturbance  or  it  nay  consist  of  a flow  of 
particles.  If  the  flux  is  in  the'  form,  of  radiant  energy,  we 
shall  not  be  concerned  with  its  frequency  distribution  o 2*  its 
plane  of  polarisation  for'  these  properties  become  important  only 
when  the  response  of  the  receiver  to  the  incident  flux  Is  considered. 
At  point  (Xq,  To,  Z0)  there  is  Indicated  an  elementary  receiving 
area  ds*  which  is  small  compared  with  all  the  dimensions  considered 
and  has  a surface  normal  H!  with  direction  cosines  ( oC  , , *")■» 

Both  the  coord. inates  of  the  elementary  receiving  area  and  the 
direction  cosines  of  its  surface  normal  are  completely  arbitrary. 

Me  shall  assume  that  each  element  of  area  of  the;  rectangular 
source  emits  the  same;  flint*  per  unit  solid  angle  in  » direction 
nciBal  to  its  surface , srd  that  it  obeys  Lambert's  Cosine  Law. 

If,  in  addition.,  we  assume'  that  the1  Intervening  space  between 
the  source  and  receiver  is  homogeneous  and  neither  absorbs  nor1 
scatters  the  flus  passing  through  It,  it  is  possible  to  calculate 

rigorously  the  total  flux  flailing  on " the  elementary  area.  da*. 

In  performing  these  calculations,  we  shall  neglect*  the  ©ontritou™ 
tions  of  multiple  reflections,  raul tipis  absorptions  and  reenissions 
between  the  surfaces  which  might  contribute  to  the  total  flux 
falling'  02  ds*.  for  ds*  will  always  be  assumed  to  be  so  snail 
or  such  an  efficient  energy  converter  that  such  effects  will  sot 
make  appreciable  contributions • 

1%.  The  radiation  'Characteristics  of  the  rectangular  source 
mast  now  be'  expressed  in  mathematical  terms.  We  postulate  that 
each  unit  area  of  the  rectangular  source  radiates  the  same  total 


♦The  units  of  flux  in  the  e.g.s.  system  are  **rgs/s«e;  In  the 
auk.  a.  system  watts. 
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**Br  total  energy  we  mean  all  the  energy  entering  within  the 
unit  solid  angle.  Thus,  if  the  flux:  Is  In  the  fers  of 
particles,  it  Bast  then  he  the  total  energy  contained  by  the 
particles  which  enter  within  this  unit  solid  angle*  If  the 
radiation  Is  periodic.  It  must  he  integrated  ©war  all 
frequencies;  .If  It  Is  transfer se,  It  unat  he  integrated 
over  all  planes  of  polarization.  The  units  of  energy  in 
the  c.goS*  system  are  ergs;  in  the  m.k.s.  system  joules. 
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energy**  B per  second  per  unit  solid  angle  in  the  direction  of 
its  normal.  On  a microscopic  tine  basis,  B,  of  coarse,  met 
fluctuate  because  for'  particles  the  flux  cannot  be  considered 
to  be  distributed  uniformly  throughout  space  and  for  waves  the 
energy  flow  must  be  periodic.  However,  if  the  times  of  observa- 
tion. are  made  large  compared  to  these  nleroscopic  disturbances, 
then  B will.  If  the  source  continues  to  radiate  In  the  sea® 
manner,  establish  an  average  -value.  In.  all  the  following  remarks 
it  will  be  assumed  that  the  times  of  observation  are  sufficiently 
large  that  such  an  average  value  fop-  B has  significance.  If  we 
consider  a snail  elementary  area  ds  of  the  rectangular  source-, 
the  total  flux  within  the  element  of  solid  angle  duo-  wist  then 
be  of  amount  B ds  da/  • We  have  suggested,  earlier  that  each 
elementary  area  of  the  rectangular  source  aast  obey  Lambert's 
Cosine  Law  (reference  (a)).  This  law  slaoiy  states  that  the 
total  flux  within  the  solid  angle1  d«J  being  radiated  by  a snail 
area  element  decreases  with  the  cosine  of’  'the  angle  between  the 
.surface  normal  and.  the  direction  of  emergence1  of  the  radiation. 
Since,  for  such  sisa.ll  areas  radiating  within  a snail  element  of 
solid  angle  du  the-  effective  cross  section  area  of  the  source 
also  decreases  with  the  cosine  of  this  sane-  angle,  B for1  a, 
surface  obeying  Lambert's  Cosine  Law  oust  be  Independent  of  the 
angle  of  emergence  of  the  flux*  A source  whose  total  flux1  falls 
off  with  the  cosine  of  the  angle  between  the  surface  normal  and 
the  direction  of  emergence-  of  the  flux  is  as  previously  stated 
known  as  a “perfect  diffuse  emitter”..  Lambert's  Cosine-  Law 
Is,  of  course,  only  an  approximation  for  many  sources  do-  not 
obey  this  law,  but  we  shall,  for  want  of  a better  description 
of  such  sources,  assume  It  to-  be  valid . The  total,  flux  radiated 
by  the  area  clement  ds  in.  a direction  © with  its  normal  within 
the  solid  angle-  du>  would,  with  the  assmsptlcm  that  it  obeys 
Lamberts  Cosine  Law,  be  expected  to  be  of  the  amount 
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15*  A similar  epaiioD  can  also  be  written  for*  the  total 
energy  emerging  in  a direction  © In  the  tire  dt  by  simply 
inserting  * dt  an.  the  end  of  the'  previous  equation.  V*  will 
not,  however,  he  concerned  with  energy  in  this  paper  far  most 
receivers  only  integrate  over  the  tine  in  the  short  interval 
in.  which  they  are  first  exposed  to  the  radiation.  After  "this 
brief  'initial  exposure  to  the  radiation,  .an  equUibrlusr  is 
established  in  which  the  total  flux  falling  on  the  receiver  is 
equol  ***  the  flux  expended  in  urn  lata  inir.g  a certain  signal 
level  plus  'the  flux  lost  by  reflection,  heat,  radiaticsi,  noise, 
and  other  such  losses  inherent  in  detectors. 


16*  If,  as  suggested  by  Pig.  12,  this  flux  is  allowed  to  be 
incident  on  the  receiving  area  ds*  at  an  arbitrary  angle  0* 
with  its  normal  K»«  we  wist,  in  order  to  determine  the  total 
flux  falling  on  da*,  limit  the  solid  angle  da  to  that  inter- 
cepter*  at  ds  by  the  area  ds**  From  the  definition  of  'the  solid 
angle,  we  can  write 


<SL* 
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a being  the  distance  between®  ana  ds®.  If  now  this  value  for 
de»  is  inserted  in  Eq.  l,  the  total  flux  which  the  elementary 
area  ds  radiates  to  the  elementary  area  ds1'  Is  then 

m 3d. s # afa'caa  #•* 
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This  equation  is  one  of  the  old  and  universally  accepted  laws 
of  photometry  (rofereru*®  (a))  and  is  commonly  referred  to  as 
halbert*  s Photometric  Law.  It  is  also  valid  for  small  sources 
which  do  not  obey  Lambert's  Cosine  Law,  providing  B is  defined 
as  the  total  energy  radiated  per  second  per  unit  area  per  unit 
solid  angle  in  the  direction  of  0* 

17*  This  equation  provides  the  starting  point  for  qpsntitatlve 
calculaticr.s  in  this  report  as  has  been  toe  case  for  others 
working  on  these  problems.  All  one  has  to  do  is  to  perform 
the  indicated  integration  over  the  surface  of  the  rectangle 
in  order  to  obtain  the  total  flux,  falling  css  'the  elementary 
receiving  area  ds*.  To  adapt  this  equation  to  the  problem  at 
hand,  it  is  necessary  to  refer  to  Fig,  12.  The  radiating 
area  eleaent  ds  can  be  expressed  in  terms  of  an  element  of 
area  of  the  rectangle,  which  obrriausly  will  be  dY  dZ.  The 
angle  0 between  the  morsel  U and  the  direction  of  a will  be 
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given  toy 


Cc*  Q 


X<^A 
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To  fix  the  magnitude  of  the  angle  9*  between  If  * and  af  we  use 
tbs  fact  that  the  cosine  of  the  angle  between  two  lines  whose 

direction  cosines  are  known  is  given  by 


cos  o'  « «*C4*S  (t,  x)  + $ Cos(eLi  /)+  /cos  (k»  z). 
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where  oC  A,  /"  are  the  direction  cosines  of  S'® 

Frcm  Figure  12  it  is  quite  a straightforward  task  to  write  down 
the  direction  cosines  of  a.  Assuming  a is  positive  when  going 


from  ds  to  ds  % we  find,  its  direction  cosines  are 
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The  magnitude  of  a is  also  easily  dieter  minted®  It  is 
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If  now  we  substitute  for  ds,  cos  9.  cos  9®  and  a in  Equation  3f 
the  integral  equation  which  must  toe  salved  in  order  to  determine 
the  total  flux  falling  on  ds  * fro®  the  rectangular  source  must 
then  toe 
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perhaps  be  noted  that  the  direction  cosines  of  the 
surface  normal  H *'  have  been  chosen  to  be  positive  in.  the  sane 
sense  .in  which  a is  positives,  Thus  H * will  always  toe  directed 
in  the  general  direction  of  the  incident  radiation,  and,  hence, 
will  toe  the  surface  normal  for  the  surface  on  which  the  radiation 
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is  not  Incidento  This  same  convention  for  Nc  has  oeen  used  for 
all  the  figures  in  this  paperQ  The  integration  over  the  Z 
direction  is  readily  carried  outQ  This  integral  is  an  integral 
of  the  form 


(a-  vijccz 

(B -Cl  +zf 


and  is  equal  to  the  sum  of  two  integrals  each  of  which  are  of 
standard  forms©  When  the  integration  was  carried  out,  the 
following  result  was  found e 
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This  equation  must  of  course  he  evaluated  at  the  limits  -w  and 
+Wo  When  this  is  dcne?  we  find  that 
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18.  To  complete  the  problem,  each  of  these  terms  must  he 
integrated  over  the  Y direction  vithin  the  limits  of  -_Z  to  +-c« 
This  looks  rather  discouraging  since  from  Rq«  8 we  seo  also  that 
A and  E are  functions  of  Y and  must  he  inserted  before  integra- 
tion is  attempted.  Some  of  the  terms  in  the  previous  equation 
will  then  he  found  to  divide  up  into  other  terms,  each  of  which 
must  he  integrated  separately.  It  is,  of  course,  impossible  to 
carry  out  such  a detailed  Integration  here,  and  about  all  that 
can  be  d CjuLCV  Is  "to  §ivs  s CmS  gonora  1 suggestions  which  might  be 
useful  to  anyone  who  cares  to  repeat  the  calculations,  and  then 
simply  give  the  final  result.  To  reduce  each  of  the  integrals 
to  a standard  form,  a change  of  variables  is  necessary.  For 
all  the  integrals,  the  same  change  of  variables  is  utilized. 

This  Is  quite  fortunate  for  all  the  integrations  can  be  carried 
out  in  terms  of  this  new  variable  and  the  various  terms  collected 
and  simplified  before  worrying  about  the  limits.  The  change  of 
variables  utilized  was 


if  » y~  Z0>  dv'  * cLZ 


(10) 


With  this  new  variable,  the  limits  on  the  integral  must  be 
changed  to  ~(Y0+/)  to  -(Y0-.Z) . In  carrying  out  these  integra- 
tions, the  following  well-known  integral  forms  were  encountered. 
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Terms  involving  arc  tangents  were  always  integrated  by  parts. 
After  these  integrations  were  carried  out,  the  terms  in  Eq.  9 
were  found  to  give  thirty-eight  separate  terms  without  inserting 
the  limits  for  v.  Fortunately,  a very  large  number  of  these 
terms  cancel  out  and  the  resultant  terras  simplify  considerably. 
The  resulting  value  for  this  integral  was  then  the  following. 
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If  now  we  insert  the  proper  limits  on  the  variable  v as sugpested 
earlier « we  obtain  the  following  equation  from  Eq.  8 for  the 
total  fiux  falling  on  the  elementary  area  ds  from  th„  rectangular 

source. 
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- * -*- amative  end  somewhat  simpler  form  of  this  equation  can 
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in  terms  of  a single  angle.  dP  may  then  also  be  written  in 
the  following  alternate  form. 
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19«  These  last  two  equations  become  quite  powerful  tools  in 
radiation  problems,  for  the  user  can  now  calculate  the  total 
flux  from  a rectangular  source  falling  on  an  elementary  receiving 
area  having  arbitrary  coordinates  (Xq^Yq, Zp)  ana  a surface  normal 
with  arbitrary  direction  cosines  tH » These  equations 

are  very  general  in  form  so  that  all  similar  equations  for  which 
the  rectangular  source  is  confined  tc  one  of  the  coordinate 
planes  at  the  center  of  the  coordinates  become  special  cases 
of  these  equations  in  which  either  the  coordinates  or  the 
direction  cosines  of  the  surface  normal  of  the  elementary 
receiving  area  have  particular  values.  After  a simple  transla- 
tion and  rotation  of  axes,  these  equations  may  also  be  used  for 
the  equally  important  case  in  which  the  elementary  receiving 
area  is  fixed  in  one  of  the  coordinate  planes  at  the  center  of 
the  coordinates,  and  the  rectangular  source  permitted  to  have 
arbitrary  coordinates  and  a surface  normal  with  arbitrary 
direction  cosines.  Using  the  equations  in  this  form,  it  is 
then  possible  to  calculate  the  total  flux  falling  on  an  elemen- 
tary receiving  area  in  the  presence  of  an  array  of  such  rectangular 
sources.  Tf  the  source  is  a surface  radiator  of  irregular  shape 
or  does  not  radiate  uniformly  over  its  surface  or  obey  Lambert's 
Cosine  Law,  it  is  usually  possible  to  divide  this  surface  up 
into  a series  of  smaller  rectangular  sources  which  radiate 
approximately  uniformly  over  their  surfaces  and  approximately 
obey  Lambert’s  Cosine  Law  within  the  limits  of  the  angles 
intercepted  at  the  source  by  the  elementary  receiving  element. 

The  radiation  fields  of  more  complex  surface  sources  can  then 
also  be  approximated  by  using  these  equations.  The  only 
precaution  necessary  in  using  these  equations  is  that  of  being 
certain  that  the  plane  containing  the  elementary  receiving 
area  does  not  intersect  the  rectangular  source,  for  if  it  does, 
these  equations  (except  in  very  special  cases)  no  longer  apply. 
Vhen  this  plane  intersects  the  rectangular  radiating  surface 
in  a line  parallel  to  one  of  the  pairs  of  sides,  corrections 
in  the  limits  can  be  made  to  take  care  of  the  reduction  in 
effective  radiating  surface  so  as  to  obtain  the  total  flux 
falling  one  one  side  of  the  elementary  receiving  element « 

Care  must  be  taken,  however,  to  shield  the  opposite  side  of 
the  elementary  receiving  area  from  radiation  so  as  to  obtain 
the  calculated  flux.  When  the  plane  containing  the  elementary 
receiving  area  intersects  the  rectangular  source  in  a line 
which  is  not  parallel  to  a pair  of  its  sides,  triangular 
sources  have  to  be  considered  and  these  were  considered  to 
be  outside  the  scope  of  the  present  report. 

20.  One  of  the  most  immediate  questions  which  confronted  us 
was  how  ve  could  be  sure  that  these  equations  are  correct  in 
view  of  the  errors  which  others  have  made  on  simpler  aspects 
o jl  uut:  j>i xn  oiuSii'  co  their  correctness,  we 

immediately  started  a search  of  the  literature  for  all  types 
of  equations  for  rectangular  sources  and  also  worked  out 
independently  a series  of  problems  in  which  (XD,  Yp,  Zp)  and 
( oC , ) had  certain  restricted  values.  This  literature 

search  was  responsible  for  the  history  of  the  subject  presented 
earlier.  The  results  obtained  by  others  and  from  these 

lb 
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independent  calculations  were  then  compared  with  the  restricted 
forms  obtained  from  Eqs.  12  or  13  when  particular  values  for 
(X0,Yq,Zq)  and  (oC  , $ , yO  were  inserted.  The  most  direct 
verification  of  the  correctness  of  these  equations  was  that 
obtained  from  Bethe's  equation.  After  a translation  of 
coordinates,  insertion  of  direction  cosines  in  place  of 
angles,  and  a change  in  sign  of  one  of  these  direction  cosines, 
all  the  various  terms  could  be  collected  and  Eq.  12  was  found 
to  be  the  result.  This  is  fortunate  for  the  solutions  of  two 
closely  related  complex  problems  involving  different  integrations 
serve  as  checks  on  each  other.  It  should  also  be  stated  that 
in  all  other  cases  in  which  independent  calculations  on  simpler 
problems  done  by  others  or  the  author  were  correct,  the  restricted 
form  of  the  general  equation  was  always  in  agreement. 

IY.  Radiation  Equations  for  Special  Simplified  Cases 


21.  Since  the  literature  on  the  rectangular  emitting  surface 
is  rather  meager,  widely  distributed  both  in  periodicals  and 
tire,  and  often  in  error,  it  may  be  of  interest  to  consider 
some  special  cases  of  the  above  formulae.  In  considering 
these  special  cases,  we  will,  insofar  as  possible,  cite  the 
originator  of  each  of  these  formulae  for  many  of  these  formulae 
have  often  been  rederived  in  the  literature.  Ve  will  also 
point  out  wherein  we  believe  errors  have  been  made  when  we 
believe  they  are  present.  In  cases  in  which  we  believe  errors 
to  be  present,  we  have  not  relied  solely  on  the  general  equation 
of  the  present  paper  to  establish  these  errors,  but  have  also 
carried  out  independent  calculations  starting  with  the  original 
problem  suggested  by  each  equation,  in  presenting  these  special 
equations,  they  will  be  considered  in  the  order  of  decreasing 
complexity. 

A.  Bethe’s  First  Formula 


22*  One  ui  the  special  cases  which  immediately  suggests  itself 
is  that  in  which  the  surface  normal  is  not  limited  in  its 
orientation,  but  the  coordinates  of  the  receiving  element  are 
limited  to  one  of  the  coordinate  planes.  Such  a formula  was 
published  by  Bethe  (reference  (n);  in  1951®  If  in  Fig.  12  we 
confine  the  receiving  element,  for  example,  to  the  i'X,Y)  plane, 
we  crust  then  in  Eq*  13  set  ZG  ^ 0.  Doing  this  we  find  that 
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This  equation  is  a very  much  condensed  version  of  Beths' s 
equation  and  will  be  found  equivalent  to  it. 


X?  Qr 


Equations  for  which  c*  , ft  , and  Y are  zero 


23 o Another  of  the  simpler  cases  in  the  order  of  decreasing 
complexity  is  that  in  which  the  elementary  receiving  area  is 
restricted  so  that  its  surface  normal  rotates  in  a plane 
parallel  to  one  of  the  coordinate  planes.  Equations  for 
these  three  cases  can  be  arrived  at  from  Eq.  13  by  simply 
setting^  - 0,  - 0,  and  / *=  0.  When  this  is  done,  the 

following  three  equations  were  found. 
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These  equations,  although  rigorous  and  apparently  new,  still, 
represent  examples  of  cases,  which  are  considered  too 
complicated  to  be  evaluated. 

Co  Equation  for  which  and  2.0  are  zero 

24.  One  of  the  simpler  variations  of  the  previous  formulae 
is  the  case  in  which  the  receiver  is  restricted  to  the  (X.Y) 
plane  in  such  a way  that  Its  surface  normal  also  rotates  In 
this  plane.  The  equation  describing  the  total  flux  falling 
on  the  receiver  in  this  case  can  be  obtained  by  setting  both 
¥ and  Z0  in  Eq . 13  equal  to  zero , or  by  simply  setting 
Zip  **  0 in  the  last  of  the  previous  equations.  In  this  case 
dF7:y  reduces  to 

r 

Jrr  Bdi  L I 
<u’Xy‘  A \Jx?rur>- 

VXe* +[/<,-* )*  XS+ us* 


^'kZyf’ty  -faui  J2iv  I 

yhS+fa+J)  * Xo*r  c/0* *£)*-  ur&\  t 


This  formula  also,  to  the  best  of  our  knowledge,  has  not  appeared 
in  the  literature. 

D„  Equation  for  which  o*.  and  Z0  are  zero 

25.  Another  interesting  case  in  that  shown  in  Fig.  13«  In  this 
case  the  receiver  is  still  confined  to  the  (X,Y)  plane  but  its 
surface  normal  -Is  restricted  to  rotation  in  the  (Y,Z)  plane. 
Clearly  from  Fig.  12,  we  would  expect  both  o<  and  ZG  to  be  zero 
for  this  caseo  This  time  we  will  substitute  these  values  in 
Eqc  12  because  the  result  is  simpler  than  the  corresponding 
result  from  Eq.  13.  Thus  we  find  that 
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Far  this  equation  to  he  valid , the  plane  containing  ds 8 must  not 
intersect  the  rectangular  sourcec  It  is  particularly  interesting 
to  note  that  when  Y0  also  Is  zero  that  dF  vanishes®  This  is  to 
he  expected  since  the  same  total  flux  woula  fall  on  both  sides  of 
the  receiver  regardless  of  its  orientation  with  the  result  that 
the  total  flux  of  a given  sign  falling  on  the  receiver  must  be 
zerCo 

E0  Jones  0 Formula  for  Figure  6 

260  Some  ox  the  simpler  formulae  in  the  literature  which  are 
closely  related  to  the  previous  formulae  are  the  Jones  formula 
for  the  case  suggested  in  Figure  6 and  the  HIgbie  formula  for 
that  in  Figure  li0  To  obtain  the  Jones  formula,  the  normal  to 
the  receiving  element  must  be  directed  along  the  radius  R drown 
from  the  center  of  the  rectangle  to  the  receiving  element® 

Clearly  then  from  Figures  6 and  12 , the  following  substitutions 
should  be  made® 


•c  * 


Y * o t 


To  obtain  the  equation  in  the  form  in  which  Jones  originally 
wrote  it,  it  is  necessary  to  substitute  in  Equation  12  rather 
than  Equations  13  and  l40  The  result  is 
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This  is  the  Jones  formula  (reference  (f))  which  was  published 
correctly  for  the  first  time  in  19110  Another  fora  of  this 
equation  can  be  obtained  by  substituting  in  Equation  14  rather 
than  Equation  12,  but  this  additions}  equation,  although  some- 
what  simpler,  does  not  have  any  particular  advantages  over  the 
form  in  which  it  was  first  published 0 

F0  Hlgbie  °s  Formula  for  Figure  11 

27c  The  Higbie  Levine  formula  for  the  case  described  by 
Figure  1.1  is  a special  case  of  dFjry  in  Equation  15c  To  reduce 
dFxy  to  the  Higbie  Levine  formula,  the  following  substitutions 
in  Figure  12  will  be  necessary,. 
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When  these  substitutions  are  made  in  dFxy  of  Equation  15*  the 
following  result  was  found  J 
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This  equation  is  in  the  form  in  which  it  was  first  published 
(reference  (i))  and  is,  of  course  correct,. 

Go  Yamanouti as  Formulae  for  Figure  8 


2b„  All  other  formulae  found  in  the  literature  were  for  cases 
in  which  the  surface  normal  to  the  receiving  element  was  parallel 
to  one  of  the  three  directions  of  the  coordinat.es0  This  is 
clearly  the  case  discussed  first  by  Yamanouti  (reference  Ch}) 
which  is  described  by  Figure  8„  To  arrive  at  the  equations  of 
Yamanouti , it  is  necessary  to  set  d-  , /d  and  jr'  equal  to  unity* 

We  cap  make  these  substitutions  in  either  Equation  12  or  Equation 
13o  •'■o  obtain  the  • qua t Ions  in  their  original  forms , it  Is 

necessary  to  use  E uation  12 0 Thus 
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These  equations  are  in  agreement  with  those  of  Yamanouti  with 
the  exception  of  a change  In  sign  of  the  w in  the  numerator  in 
the  coefficient  of  the  second  term  in  dFv,  and  a change  in^si^n 
in  the  .£  under  the  ■nodical  in  the  dsnomximlGx'  wf  the  uweiflolcivt 
of  the  second  term  in  dF  0 Both  of  these,  we  believe,  are  mis=” 
prints  and  we  have  written  to  Yareanouti  concerning  this0 


Ho  Condensed  Form  of  Yamanouti Js  Equations 

29o  It  should  also  be  mentioned  that  a condensed  version  of 
these  equations  can  be  obtained  by  making  substitutions  in 
Equation  13*  These  equations  are 
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I0  Equations  of  Moon  and  Spencer 

1 

30o  The  closest  equations  to  those  of  Yamancuti  are  those  of 
Moon  and  Spencer  (reference  Sm))0  In  fact,  the  problem  discussed 
by  Moon  and  Spencer  is  identical  with  that  of  Yamanouti,  but  the 
manner  in  which  the  equations  she  sot  up  is  a little  different® 
Refening  to  Figure  14,  Moon  and  Spencer  determined  the  coordinates 
of  the  receiving  element  from  one  of  the  corners  of  the  rectangle. 

To  transform  Yamanouti 9s  equations  to  those  of  Moon  and  Spencer, 
a translation  of  coordinates  is  necessary®  Referring  to  Figure  14, 
the  translation  must  be  such  that 
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Introducing  this  translation  of  coordinates  into  Yamanouti “s 
equations,  we  find  that 


dFx  - B&lL  -ac-  L-,  Ur*,  , , . 


^ rrr 


+ U*-' , 


i 'JKFRZstv)* 


V -=-  , 

\ZXo*r*£  } fx7F~s.  " ~rm  P== 

i *x°  *-■“■ 


: ' — / 

^^ZrirP'J 


XoAt^ 


3 -X  f 

iWHiir*)*)  7?w^, 


Jyit,}  -/  d *r  _ —£?C  — | 


V A"./-  foX-S  ^ f"  (i  o JL,  A*"  _ \ 

' j J%  j f^TSTT  2>  y / 


U~,~j  3ttr  i~AT  . ,-iT /'{ 

/x/^&TiZ)^  0 fyc^2~i+^  ( / v 


kM^7^  ( / 

yf 


; 9 

f.v^  — ’f  <**■  7^ 


JZFl-nr* 


=£=====■  j / . , 


. ,vk 


NAVOJRD  Report  29&0 


These  equations  are  the  equations  of  Moon  anti  Spencer  (reference 
(m))  with  the  exception  of  the  B term  which  they  define  in  a 
different  manner „ 

J,  HigbieBs  Formula  for  Figure  9 

31*  An  equation  which  is  equivalent  to  one  of  the  Yaraanouti 
equations  is  that  derived  by  Higbie  (reference  (i))  suggested 
in  Figure  9»  To  transform  Equation  12  to  the  case  suggested  by 
Figure  9?  the  following  substitutions  are  necessaryQ 
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Making  these  substitutions  in  equation  12 9 


This  is  Higbie 9s  equation  published  in  1925  and  it  is?  of 
course,,  correct „ This  equation  is  equivalent  to  dFy  in 
Yamanouti 5s  equatlonso 
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K»  Bethe*s  Formula  for  Figure  15 

32 o Another  recent  variation  of  Yamanouti’s  equations  Is 
that  given  by  Be the  (reference  (n))  shown  schematically  In 
Figure  15 » In  this  case,  the  surface  of  the  receiving  element 
is  parallel  to  the  rectangular  emitting  source  and  Its  motion 
Is  confined  to  a plane  determined  by  the  Xc  axis  and  the 
line  Vo  — Z0  * This  is  a special  case  of  the  first  of 
Yamanoutl' s equations  (that  is,  withc<  = / ) in  which  y0  Is 
simply  set  equal  to  c In  order  to  write  this  equation  in.  the 
final  form  in  which  it  appears,  it  is  preferable  to  make  these 
substitutions  directly  in  Equation  13  o Doing  this.,,  dFx  has  the 
following  value o 
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This  is  Be the’ s equation  in  the  form  in  which  it  was  written. 
It  is  more  useful  in  connection  with  square  sources  since  in 
this  case  the  flux  dFx  would  he  that  hi  a plane  normal  to  the 
square  surface  and  through  its  diagonal. 


1 


3r 


NAVOED  Report  29&0 


Lo  Jones  ’ Formula  for  Figure  5 

33c  One  of  the  simpler  special  cases  of  the  Yamanouti 
equations  is  that  suggested  by  Figure  5 in  connection  with 
Jones  3 contributions o In  this  case?  the  receiving  element  is 
confined  to  a plane  parallel  to  the  C X,  y ) plane  but  displaced 
a distance Z0  ~ and  the  surface  normal  lies  in  a direction 

parallel  to  the  y axis0  Obviously  we  must  then  make  the 
following  subs ti tut ions o 


Z0=  . 


Making  these  substitutions  in  dPy  of  Equations  20,  we  *'ind  that 


cLFy 


r 


X* 


AV 


(^.4-  — - / ^ — Cz&~—  _ A& ~f-n.Tr1  *v  | 

(25)  •* 


Except  for  the  constant  factor  outside  which  was  neglected  but 
could  be  inferred  from  previous  equations  and  the  misprint  of 
an  exponent  as  a subscript,  this  is  the  Jones"  formula  (reference 
Cf ) ) .•  It  has  often  been  rederived  in  the  literature,. 

Mo  Higbie “s  Formula  for  Figure  10 

3*+,,  Another  special  case  of  these  general  equations  is  that 
of  Higbie  and  Levine  represented  schematically  by  Figure  10o  In 
this  case,  the  receiving  element  is  located  at  a distance  a in 
their  direction  from  a normal  through  one  corner  of  the  rectangle 
and  is  oriented  with  its  surface  parallel  to  that  ©f  the  rectangle0 
In  order  then  to  transform  these  general  equations  t©  Higbie  “s 
formula,  the  following  substitutions  must  be  made0 


^ i X>  ® ~£j  Z0  - CL  -Oj~,  £/  e Z ZuT  ^.  sTr-x-'  • 


The  most  direct  method  of  arriving  at  Higbie “s  equation  is  that 
of  using  Equation  13<=  Making  these  substitutions  in  this  equation,. 
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(26) 


This  is  the  Higbie  Levine  formula  (reference  (i))  vhich  appeared 
first  in  1926  and  it  is,  of  course,  correct., 

N.  Jones'  Formulae  for  Fig*  3 and  Fig*  4 

jp»  Jones  aj.sc  discussed  -cwo  aQuicicuax  cases  wnj.cn  are  repre- 
sented by  Figs.  3 and  4.  Obviously,  if  we  are  to  transform  the 
general  case  in  Fig«  12  to  that  suggested  by  Fig.  3j  the  following 
substitutions  must  be  made* 


<=£  ~ >,  y0  - 'st,  z. 


mu  t xi  r. — x s v j>  . M *1  O *1  •» 

incat:  auuatxiiiAiiuiih  v«ii  uc  .tiauc  jk.it.  cua^t  kjh\s  wa  ^400  ac j aj? 


20,  21  with  the  same  result, 

{tFy  = 

i ££ 

y/X*rf&Jl)x 


sj 


7~a-r ' 


^5  ^ 
^-4 


sFun  ~f 


&t*r~ 


Jim, 

JxTZlii!)* 


(27) 


This  equation  was  first  published  (reference  (f))  incorrectly 
is  1909  but  was  corrected  in  a following  paper  in  1910.  For 
Fig.  4,  the  receiving  element  is  parallel  to  the  rectangular 
source  and  is  located  on  a surface  normal  through  the  center 
of  the  rectangle.  This  again  is  a special  case  of  the  previously 
mentioned  equations  in  this  paragraph  in  which  K 1 and  and 
are  both  zero.  Making  these  substitutions, 
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df~}e  38  £8<U*  I ~ur- 


A. 


w?***  vfcW*  v/f2^£ 


7^X77 


|(2b) 


This  equation  is  again  correct  and  appears  in  Jones5  paper 
(reference  (f)}  in  1911 0 

0C  Bordoni 8s  Formulae  For  Figure  1 and  Figure  2 

360  The  last  two  equation?  which  we  shall  mention  are  those  of 
Bordoni  (reference  (e>)  represented  by  Figures  1 and  2C  These  are 
slight  modifications  of  the  cases  already  discussed  in  on^or.Mwi 
with  Jones5  Danersr,  Thu?  to  reduue  our  general  equations  to  those 
represented  in  Figure  2,  the  following  substitutions  must  be  made 
in  Equation  120 


*-4  Z. 


O 


Doing  this. 


A 


dfy  = /Sds'jh*-!  £ 


tm-t 


^<r 


A/.Z  ./a 
ys\#  rx^c-Jr/ 


C’29) 


To  obtain  the  aquation  representing  Figure  1,  the  following 
substitutions  must  be  made0 


~ yi  ■■  2-0 


5 ' 


V 


Thus  Bordoni 8s  equation  for  the  case  in  Figure  1 is 


— /t^s  /\  ■ — -bm'f  £ 


/XfrZr** 


o2j£ 


=±v^.  ?' 
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As  suggested  earlier,  both  of  these  equations  are  correct  and 
apparently  represent  the  first  contribution  to  the  radiation 
fields  of  rectangular  sources  of  finite  dimensions* 

P0  Special  Photometric  Equations 


37*  There  are  two  other  well-known  simple  photometric  lav/s 
which  can  be  used  al30  as  a check  on  these  equations.  Since 
Lambert’s  Photometric  Equation  given  earlier  as  Lq,  3 17ns  used 

it  also  ought  to  be  contained  in 


in  deriving  these  equations,  i 
these  equations  as  a special  c 


rp  Oi  1 4 M l-t  r»  4*  4 n 4 M 


the  case,  it  is  necessary  simply  to  assume  that  Xo  is  much 
larger  than  IQi  Z0,  £ , and  w.  If  ue  make  this  assumption  in 
Eq*  13  and  expand  the  arctangent  terms  in  series  retaining 
only  the  first  terms,  we  find  that 


C*f  £ 


K 0 A 


(3D 


which  is  in  accordance  with  Lambert’s  Photometric  Law-  Finally 
if  we  tnake^  and  v very  large  and  of  the  same  order  of  magnitude, 
then  for  increasingly  smaller  values  of  XOJ  Y0s  and  Z0,  Eq0  13 
reduces  to 


~ ir  Sets  ' ^ (32) 

again  a well-known  result* 

36  * It  should  perhaps  be  mentioned  again  that  caution  must  be 
exercised  in  the  use  cf  these  equations.  If  the  plane  con- 
taining the  elementary  receiving  area  intersects  the  rectangular 
source,  those  equations  are  then  no  longer  applicable  except 
in  very  limited  cases 0 

V,  Radiation  Field  of  a Square  Surface  Source 
A,  Deviations  from  Inverse  Square  Law 

39*  As  a particular  example  of  these  equations,  a few  calcu- 
lations were  carried  out  for  a square  surface  source.  One  of 
the  simplest  and  most  often  calculated  cases  for  these  problems 
is  that  in  which  the  elementary  receiving  area  is  placed  on  tha 
X0  axis  parolloi  to  the  source  ( c<  =- 4 Xo  ■*  Za  o ),  For  such 
a case,  our  general  equation  goos  over  into  the  following 
simple  form: 


c£/~jc  — f 


L 


! v A<r 

29 


7^0 


sJXj 


1 


(33) 


NAVOBD  Report  2980 

v being  the  width  of  the  square*.  In  Fig.  16  we  have  plotted 
some  data  arrived  at  from  this  equation.  In  this  figure,  the 
ordinate  is  dF«/B  dg‘  and  the  abscissa  is  XD  in  units  of  the 
width  w.  The  full  line  curve  in  this  figure  gives  the  depend- 
ence of  dFx/B  ds*  on  X0  as  would  be  expected  from  this  equation. 

40.  Often  in  problems  of  this  sort,  the  inverse  square  law  is 
applied  without  being  sufficiently  careful  as  to  the  magnitude 
of  the  error  introduced  by  its  use  . To  indicate  the  type  of 
agreement  to  be  found  between  the  exact  Eq.  33  for  the  square 
and  that  tc  be  expected  on  the  basis  of  an  inverse  square  law, 
we  ha%'e  included  a second  dotted  curve  in  Fig.  16  showing  the 
values  of  dFx/B  els’  found  from  the  inverse  square  law.  If  we 
set  aC  = l in  Eq.  31,  the  equation  describing  the  dependence 
of  dFx  on  X0,  assuming  an  inverse  square  law,  must  then  be 

OLzlH 

X V*  A * 

r • o 

For  small  values  of  XOJ  the  inverse  square  law  gives  larger 
values  than  the  exact  equation,  the  divergence  between  the  two 
values  being  larger  as  Xq  decreases.  For  small  values  of  Xq 
compared  with  w,  the  ordinate  for  the  full  line  curve  approaches 
■JT  while  that  for  the  dotted  curve  gees  to  infinity  v 

Bo  Special  Case  for  which  - Y0  « Zc  = 0 

41.  Another  quite  simple  example  is  that  for  which  the  elementary 
receiving  area  5s  still  on  the  Xq  axis,  but  its  surface  normal 

Is  free  to  rotate  in  the  (X,  Y ) plane.  For  such  a case  in  our 
general  equations,  we  must  set  y'  = 0 and  Y0  * ZD  * 0.  Doing 
this  we  find  that 

dfm  Sds  ' )'  An-/  =—==■]  (3?) 

( fXf+ur * , 

This  equation  is  identical  in  form  with  E-q.  33  except  for  the 
presence  of  the  ^ „ If  we  plot  dF/B  ds*  as  a function  of  Xq 
and  we  can  then  represent  these  values  by  polar  plots  or 
the  form  shown  in  rig.  iy.  For  a fixed  value  of  XD,  dF/B  ds’ 
falls  or.t'  simply  with  the  cosine  of  the  anglo  9' . “in  this 
figure,  the  lengths  of  the  vectors  at  each  of  the  indicated 
values  Xc  (X0  35  w,  2w,  3w)  contained  in  the  circles  of  diameters 


V r(*-r/ 
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are  the  magnitudes  of  dF/B  ds°,  K having  the  values  I,  2,  and  3 
for  the  Indicated  values  of  Xge  Below  each  of  the  circles  is 
indicated  the  maximum  value  of  dF/B  ds*.  In  each  case,  th© 
elementary  receiving  area  itself  is_normal  to  these  vectors. 

For  each  value  of  Xos  values  of  dF/JB  ds‘  are  significant  only- 

up  to  the  angle  ©'  at  which  the  plane  of  the  elementary  receiving 

surface  Itself  begins  to  intersect  the  square  surface 

(for  X0  = v,  ©9  *’45°;  X0  * 2v5  ©*  = 60°,.  Beyond  this  angle, 

these  calculations  of  course  have  no  value* 


C.  Isophotopic  Plots  for  (Y,  Z)  Plane  with  ©C  * 1 


42,  In  the  last  series  of  figures,  we  have  plotted  field  data 
of  a mere  detailed  type.  For  each  of  these  figures,  the 
elementary  receiving  area  was  confined  to  a given  (Y,  Z)  plane 
(that  is,  having  a fixed  XQ  value)  and  had  a surface  normal 
parallel  to  the  X0  axis  (that  is  cA.  * 1),  In  carrying  out 
these  calculations,  polar  coordinates  were  introduced  in  the 
(X,  plane  and  calcinations  were  made  for  a rixed  angle  (jf  with 
the  ZG  axis  and  at  various  distances  from  the  Xq  axis.  If  in 
the(Y.  Z)  plane  the  radius  r and  the  distance  Xc  are  expressed 
in  units  of  w,  Eq.  24  can,  for  purposes  of  calculation,  fca 
written  in  the  form 


yy-__  6ds*  r to*  ef-/ 

* A JXa  z+(r 
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In  arriving  at  numerical  values  from  this  equation,  values  for 
dFx/B  ds'  were  calculated  along  fixed  directions  (that  is, 
constant  and  variable  r)  in  each*  of  six  (Y,  Z)  planes  for 
which  X0  had  the  values  1,  2,  4,  6,  8,  and  10.  Intervals  of 
10°  in  <9  were  used  in  each  of  these  planes.  From  this  data, 
curves  were  plotted  in  eac-h  plane  along  each  of_the  chosen 
directions  (that  is,  for  fixed  4>  and  variable  r)  and  values 
of  r chosen  for  which  the  quantity  dFx/B  &s*  had  certain 
constant  values.  These  data  were  then  plotted  on  sheets  of 
coordinate  paper  corresponding  to  each  of  these  six  planes, 
and  all  the  points  along  each  of  the  directions  for  which 
dFx/B  ds*  had  the  same  value  were  joined  together.  The  results 
of  these  detailed  calculations  and  plotting  are  giver*  in  Figs. 

18,  19,  20,  21,  22,  and  23.  In  each  of  these  figures,  the 
values  of  dFx/B  ds'  along  each  of  these  isophot&pic  lines  is 
indicated  directly  on  the  figure.  The  square  itself  is  plotted 
to  scale  at  the  center  of  each  of  these  figures.  Most  of  these 
isophotopie  lines,  except  for  Fie.  18.  can  be  approximated  quite 
wen  Dy  circxes  ix  tho  calculated  and  plotted  values  are  accurate 
to  three  figures. 

43.  Figures  of  this  type  can  be  quite  useful  in  making  rapid 
determinations  of  the  total  flux  or  flux  density  on  an  elementary 
receiving  area  placed  in  various  planes  and  having  surface 
normals  with  various  orientations.  All  one  needs  to  know  are 
the  coordinates  of  the  point  In  question  (which  are  determined 
from  the  figure  directly),  the  magnitudes  of  B and  ds*  and  then 
either  dF  or  dF/d3’  can  be  determined  directly  from  such  curves. 

VI.  CONCLUSIONS  AND  ACKNOWLEDGMENTS 

44.  In  this  report,  we  have  presented  a rather  detailed 
investigation  of  the  subject  of  the  radiation  characteristics 
of  rectangular  sources,  when  these  sources  radiate  uniformly 
over  their  surfaces  and  obey  Lambert's  Cosine  Law.  This  has 
been  done  from  a very  general  viewpoint  so  that  the  elementary 
receiving  area  may  have  arbitrary  coordinates  end  a surface 
normal  with  arbitrary  direction  cosines.  After  a simple 
translation  and  rotation  of  axes,  these  equations  may  also  be 
used  for  the  equally  important  case  in  which  the  elementary 
receiving  area  is  confined  to  one  of  the  coordinate  planes  at 
the  canter  of  the  coordinates  and  the  rectangular  source 
permitted  to  have  arbitrary  coordinates  and  a surface  normal 
with  arbitrary  direction  cosines.  Using  the  equations  in  this 
form,  it  is  then  possible  to  calculate  the  total  flux  failing 
on  an  elementary  receiving  area  in  the  presence  of  an  array  of 
rectangular  sources  of  various  sizes  having  arbitrary  coordinates 
and  surface  normals  with  arbitrary  direction  cosines,  for  the 
total  flux  will  simply  be  equal  the  sum  of  the  individual 
contributions  of  each  of  these  sources. 
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45®  These  equations  are  also  directly  applicable  to  problems 
for  wbieh  the  previously  designated  elementary  receiving  area 
becomes  the  source  of  the  flux  and  the  rectangular  surface  the 
receiving  surface.  Thus  one  can,  by  utilizing  these  equations, 
calculate  the  total  flux  entering  a “rectangular  aperture  from 
a small  elementary  source  area  when  the  source  or  the  receiving 
area  has  arbitrary  coordinates  and  a surface  normal  with 
arbitrary  direction  cosines®  Assuming  that  the  source  in  each 
case  radiates  uniformly  over  its  surface  and  obeys  Lambert’s 
Cosine  Lav,  the  only  limitation  in  the. use  of  these  equations 
is  that  of  being  certain  in  each  case  that  the  plane  containing 
the  elementary  source  or  receiving  area  does  not  intersect  the 
rectangle,  for  if  it  does,  these  equations  (except  in  very 
special  cases)  are  no  longer  applicable. 

46.  A search  of  the  literature  has  aisoheen  made  which 
resulted  in  the  finding  of  seventeen  formulae  describing 
various  aspects  of  the  radiation  fields  of  rectangular  sources. 
V.c  have  pointed  out  th«t  all  these  formulae  are  apparently 
very  different,  and  that  there  are  a few  errors  both  in 
transcribing  and  arriving  at  these  formulae,,  We  have  shewn, 
however,  that  each  of  these  formulae  (except  Be the’ s)  can  be 
arrived  at  from  the  more  general  formulae  presented  in  this 
report  by  simply  substituting  particular  values  for  X0»  Yo, 

ZQ,  e*4-  » & » and  o * ^Tien  discrepancies  appeared  between 
the  formula  arrived  at  from  the  general  equation  presented  in 
this  report  ana  the  corresponding  equation  published  by  others, 
we  have  shown  where  this  discrepancy  exists  and  have  to  our 
satisfaction  established  these  discrepancies  as  errors  in  the 
equations  of  others.  As  a result  of  this  investigation,  we 
believe  that  all  of  these  formulae,  with  the  exception  of  Jones’ 
formula  for  Fig.  5 and  the  formulae  of  Yamsnouti  for  Fig.  8. 
are  correct  when  they  are  written  in  the  final  form  in  which 
the  authors  intended.  The  errors  appearing  in  Jones'  formula 
are  readily  found  from  preceding  equations,  but  the  errors 
in  signs  made  in  Yaraanouti's  equations  cannot  be  inferred 
from  his  paper  because  they  involve  complex  integrations  which 
cannot,  because  of  their  length,  be  discussed  in  detail. 

These  errors  in  sign  may,  however,  be  the  result  of  a misprint 
or  due  to  an  error  in  copying.  It  perhaps  should  also  be 
mentioned  again  that  the  Jones'  formula  for  Fig.  3 appeared 
incorrectly  in  his  first  publication  but  was  corrected  in  a 
following  publication.  Waldram's  formula  was  the  only  one 
which  was  not  discussed,  for-  his  formula,  by  his  own  admission, 
was  only  an  approximate  one. 

V?.  The  equations  presented  in  this  report  for  the  rectangular 
source  appear  to  be  new.  The  first  form  of  the  equation, 
although  not  specifically  written  by  Yamanouti,  can  nevertheless 
be  obtained  from  Yamancuti's  components  of  the  "light  vector" 
if  the  light  vector  is  written  for  an  arbitrary  position  and 
orientation  and  certain  definitions  are  brought  into  agreement. 
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As  mentioned  before,  however,  a discrepancy  in  two  signs  appears 
betv/een  the  two  equotionsc  The  first  form  of  the  equation 
presented  in  this  report  will  also  be  found  to  agree  with  Bathe 8s 
much  longer  equation  after  a translation  of  coordinates  insertion 
of  direction  cosines,  change  in  sign  of  one  of  these  direction 
cosines,  and  a collection  of  the  various  termso  The  alternate 
form  of  the  equation  seems,  however,  t^  be  new, 

4b o In  the  last  section  of  the  paper  we  have  presented  a few 
simple  examples  of  calculations  for  the  radiation  field  of  a 
square  surface  source  and  also  some  more  detailed  plots  of  the 
radiations  field  when  the  elementary  receiving  area  is  parallel 
to  the  square  source  and  is  confined  to  planes  normal  to  its 
axis*  Calculations  of  the  latter  type  are  very  long  and  require 
mechanized  equipment  to  carry  them  outc  Curves  of  this  type  are, 
however,  very  useful  for  rapid  evaluations  of  the  total  flux  or 
the  flux  density  on  the  elementary  receiving  area0  The  number 
of  isophotopic  curves  which  one  plots  would  depend  on  the 
accuracy  and  variability  in  the  position  of  the  elementary 
receiving  area  which  could  only  be  anticipated  by  the  user, 

49r  The  examples  which  ve  have  cited  do  not  illustrate  the 
full  versatility  of  these  equations,  for  they  can  be  used  to 
determine  the  total  flux  or  the  flux  density  on  an  elementary 
receiving  area  which  is  restricted  to  any  type  of  surface. 

All  one  needs  to  solve  the  particular  problem  are  the  coordinates 
Xo9  Yos  ZQ  and  the  direction  cosines  of  the  surface  normal  over 
tne  complex  surface,  and  then  isophotopic  lines  can  be  plotted 
over  this  surface  in  much  the  same  way  as  was  done  in  Figs*  lb 
to  23,  Radiation  fields  of  targets  having  complex  geometries 
can  also,  with  sufficient  patience,  be  calculated  with  these 
equations,  for  there  are  many  targets  which  can  be  approximated 
by  a series  of  rectangular  sources.  Mechanized  computing 
machines  will  be  quite  helpful  in  carrying  out  these  calculations, 
for.  without  them,  these  equations  become  quite  tedious  to 
evaluate  even  for 'single  isolated  points, 

50o  This  report  may  be  subject  to  criticism  because  it  is 
completely  theoretical  and  has  no  experimental  data  to  justify 
these  equations.  Further  objections  may  also  be  raised  because 
some  surfaces  may  not  radiate  uniformly  over  their  surfaces 
and  may  not  obey  Lambert es  Cosine  Law,  To  the  first  objection 
we  can  only  say  that  detailed  experimental  measurements  on  such 
problems  require  a considerable  amount  of  time  aid  this  time 
simply  was  not  available  Even  after  measurements  are  taken, 
they  may  not,  because  of  the  second  objection,  describe 
adequately  the  theoretical  calculations0  In  cases  of  this 
sort  the  measurements  are  further  complicated  because  one  must 
also  determine  how  B depends  on  the  coordinates  ana  on  the  angle 
of  emergence*  Knowing  the  dependence  of  B on  the  coordinates 
and  on  the  angle  of  emergence-  additional  integrations  must  be 
attempted  to  solve  the  problem  rigorously*  If  these  cannot 
be  done,  the  radiating  surface  must  then  be  broken  up  into  a 
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series  of  smaller  surfaces  for  which  these  conditions  of 
uniformity  and  satisfying  of  Lambert’s  Cosine  Law  are  valid. 

51 o The  author  is  indebted  to  Gerald  Fine  of  this  Laboratory 
who  expressed  a considerable  interest  in  this  problem  and 
carried  out  detailed  calculations  which  were  necessary  for 
plotting  Figs.  18,  19,  20,  21,  22,  and  23.  These  calculations 
- were  very  long  ana  tedious  and  wore  done  with  the  mechanized 

computing  facilities  of  the  Naval  Ordnance  Laboratory. 

4 Harold  Feldman  of  this  Laboratory  plotted  these  data  and  is 

responsible  for  the  drawing  of  these  figures.  The  author 
would  also  like  to  acknowledge  his  interest  in  this  work. 
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FIG.  I 

BORDONl’S  FIRST  PROBLEM 


FIG.  2 

BORDONi’S  SECOND  PROBLEM 
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FIG.  3 

JONES'  FIRST  PROBLEM 
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JONES1  SECOND  PROBLEM 
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FIG.  5 

JONES'  THIRD  PROBLEM 
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JONES'  FOURTH  PROBLEM 
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FIG.  7 

WALDRAM'S  PROBLEM 


FIG.  8 

YAMANOUTi'S  COMPONENTS 
OF  THE  "LIGHT  VECTOR" 
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HIGBiE'S  FIRST  PROBLEM 


FIG.  10 

HIG8IE'S  second  problem 


FIG.  M 

HIGBIE'S  THIRD  PROBLEM 


FIG.  12 

GENERAL  PROBLEM  SOLVED 
IN  THIS  REPORT 
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ELEMENTARY  RECEIVING  AREA  CONFINED 
TO  (X,Y)  PLANE  WITH  SURFACE  NORMAL, 
FREE  TO  ROTATE  IN  (Y,Z)  PLANE 


FIG. 14 

SPENCER'S  AND  MOON'S  MODIFICATION 
OF  THE  GENERAL  PROBLEM 
PRESENTED  BY  YAMANOUTI 
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F5G.  15 

ONE  OF  BETHE'S  SPECIAL  PROBLEMS 


FIG.  17 

ELEMENTARY  RECEIVING  AREA  FOR  GIVEN  VAL 
SURFACE  NORMAL  IS  ROTATED  IN  (X,Y)  PLANE 
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FIG.  21 

radiation  field  of  a square  source 

FOR  = ! AND  Xo  = 6w 
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